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WAVE- INDUCE)  MOTIONS  OP  A  LARGE  ROCKET  VBfICLE  DRIFTING  IN  A  VERTICAL  ATTITUDE 


BY  J.  J.  Leendertse 


ABSTRACT 

Pesults  are  presented  of  a  theoretical  and  experimental  laboratory 
study  of  the  movements  of  a  large  solid -propellant  rockat  vehicle  drift¬ 
ing  in  a  vertical  attitude  in  uniform  waves  and  in  the  wave  environment 
of  the  open  sea.  Measurements  were  made  of  the  movement  in  heave  and 
pitch  of  a  1-120  s 'ale  model  of  a  17‘>ft  long  vehicle.  Experimentally  ob¬ 
tained  results  compared  well  with  thoae  obtained  analytically. 

It  appears  that  heave  and  pitch  are  linear  functions  of  the  wave 
height  and  nonlinear  functions  of  the  wave  frequency.  Pitch  and  surge 
are  coupled. 

On  basis  of  the  linear  response  of  tne  vehicle  to  uniform  waves,  sta¬ 
tistical  features  of  the  response  in  neave  and  pitch  to  ocean  wavea  with 
a  particular  spectrur  are  calculated. 
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I.  IHTRODUCTION 


The  object  of  this  study  was  to  determine  tne  motions  of  a  large 
rocket  vehicle  drifting  vertically  in  the  vave  environment  of  the  open 
sea  prior  to  a  sea  launch.  The  waves  in  this  environment  (hereinafter 
called  the  "seaway")  vary  considerably  in  height,  shape,  and  period,  de¬ 
pending  on  the  wind,  the  water  depth,  and  the  wind-al x’ected  s*n  «c* 
fetch. 

At  present  tne  best  representation  of  the  seaway  can  be  made  by  use 
of  a  mathematical  model  of  a  random  process,  which  presents  no  values  of 
the  environment  as  a  function  of  time,  but  presents  statistical  values 
euch  as  averages  and  the  frequency  of  occurrence  of  a  certain  water  level. 
In  this  approach  the  seaway  is  taken  as  the  sum  of  a  large  nuadber  of  uni¬ 
form  wave  trains,  each  with  different  amplitudes,  and  the  profiles  of  the 

individual  vave  trains  are  assumed  to  be  sine  curves  according  to  Airy's 

(5) 

tneory. v  ' 

Ihe  magnitude  of  toe  amplitudes  of  tne  vave  traine  and  their  distri¬ 
bution  over  a  range  of  frequencies  can  be  calculated  by  a  standard  method 

(3) 

from  the  wind  velocity,  wind  duration,  and  the  fetch.  '  The  result  of 
such  a  calculation  can  be  presented  In  the  fora  of  a  wave  spectrum,  which 
is  a  distribution  of  the  mesa  squares  of  vave  amplitudes  of  the  seaway  la 
a  given  increment  of  the  frequency  over  the  vave  frequencies.  Wave  spectra 
in  the  open  sea  can  alf.o  be  obtained  directly  from  records  of  floating  or 

,  .  •  (6,7,9) 

snlpbcrne  ware  recorder:;. 

The  study  of  trie  motions  of  the  vehicle  in  a  seaway  described  by  a 
particular  spectrum  van  mode  in  two  parte.  la  the  first  part  the  responses 
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of  the  vehicle  to  uniform  waves  with  different  heights  and  frequencies  were 
studied  theoretically  and  experimentally.  For  the  experimental  work,  the 
movements  of  1-120  scale  models  of  the  vehicle  were  measured  simultaneously 
vith  the  wave  in  a  wave  tank  (Fig.  l). 

In  the  second  part  of  the  study,  the  responses  of  the  vehicle  to  sea¬ 
ways,  the  latter  represented  by  mathematical  models  of  the  random  process 
of  the  ocean  surface  elevation,  w-ro  stu^iM,  and  statistical  features  of 
the  responses  to  a  particular  seaway  were  obtained. 


II.  THEORY 
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THEORY  OF  MOVEMENTS  INDUCED  BY  UNIFORM  WAVES 
General 

The  movements  oi'  a  rocket  vehicle  drifting  in  a  wave  environment  can 
be  described  by  considering  the  movements  of  its  mass  center  relative  to 
a  rectangular  coordinate  system.  In  analogy  to  the  movements  of  a  ship 
in  a  seaway,  the  movements  of  the  vehicle  in  its  six  degrees  of  freedom 
are  taken  as  foil ova  (Fig.  2)  : 


heave  ■=»  vertical,  motion  of  the  mass  center 


surge  =<  horizontal  motion  of  the  mass  center  In  the  direction  of 
the  wave  propagation  (drift) 

evay  -  horizontal  motion  of  the  mass  center  normal  to  the  direc¬ 
tion  of  the  wave  propagation 


pitch  *  angular  motion  in  a  vertical  plane  through  the  direction 
of  the  vave  propagation 

roll  =•  angular  motion  in  a  vertical  plane  normal  to  the  direction 
of  the  vave  propagation 


ynv  -  angular  motion  around  the  vertical  axis 


Since  the  vehicle  studied  is  symetrical  in  any  horizontal  croso  sec¬ 
tion,  it  may  be  assumed  teat  coupling  exist*)  only  between  the  vave  train 
anu  the  modes  of  movement  in  the  vortical  plane  through  the  direction  of 
•-rave  propagat i .-a :  namely,  heave,  surge,  and  pitch. 

For  the  purpoee  of  analysis  tne  vehicle  is  approximated  by  a  aeries 
of  cylinders  of  different  dimensions,  as  shown  in  rig.  3#  and  is  assoaec 
to  be  drifting  with  zero  speeo  in  a  uniform  wave  train. 


'  -11  •' 
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The  wave  profile  is  approximated  by  a  sine  curve,  and  tne  water  par¬ 
ticles  are  assumed  to  describe  circular  or  elliptical  orbits,  depending 
on  the  wave  length-depth  ratios.^ 

The  origin  of  the  coordinate  sy6t«n  that  will  be  used  to  describe  the 
waves  and  the  movements  of  the  vehicle  is  taken  at  the  still-water  level 
in  the  axis  of  the  series  of  cylinders .  The  I  exit,  iu  taken  positive  In 
cue  direction  of  the  vare  propagation,  and  the  Z  arie  is  vertical,  sad 
positive  upward. 

The  horizontal  and  vertical  displacements  of  a  particle  from  it6  mean 
position  located  at  a  distance  2  under  the  still-water  surface  are 


A 


coeh  k(d+i) 
sinh  kd 


sin  (kx  -  cot  -  c) 


A  -  <*  - «) 


where 

*  horir.ouisl  displacement  of  a  particle 

u 

•  ^  -  vertical  displacement  of  a  particle 

A  --  amplitude  of  the  va re,  which  is  equal  to  half  the 
wave  height  K 

If.  a  2"  i\ 

\  *  wave  length 

c  *-  circular  frequency  (=  2r  7) 

c  -•  phase  angle 

t  -*  time 

1  r  vave  period 

d  water  depth 


(1) 

(£) 


ilea  vo 
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Considered  initially  in  this  analysis  are  tne  movements  of  one  of  the 
cylinder*  of  the  vehicle-  It  is  assumed  that  the  presence  of  the  sec¬ 
tions  of  the  vehicle  in  the  fluids  doos  not  change  the  original  pressure 
distribution  in  the  wave,  according  to  the  so-called  Fronde -Krylov  hypo¬ 


thesis 


(10) 
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If  Z  is  the  upward  displacement  of  the  mass  center  of  the  section, 

the  equation  of  not'*  ;n  can  be  written  ir.  the  manner  presented  by  Wilson 

for  the  movement  of  ships  in  regular  waves  as 

F  (p)  +  S  -  HZ  ♦  M"  ( Z  -  W  )  t  1/2  „  C  F  (W  -  Z)2  .  M  g  (5) 
z  n  i  n  d  n  n  n 

n  n 

where,  for  the  particular  section 


h  <p> 

n 


v" 

z 


total  vertical  upward  force  froct  the  wave  pressures 
the  horizontal  planes  of  the  section 

total  connection  force  with  other  cylinders 

mass  of  the  section 

added  mass  in  heaving  motion  (nth  plane) 


drag  coefficient 


?  •  significa-  c  cross  section  of  the  section 


density  of  water 


g  -  acceleration  due  to  gravity 

■  -  average  of  the  vertical  velocity  cf  tr.e  water  on  the  nth  ,lar.e 


W  -  average  cf  the  v*«ticai  acceleration  cf  tne  water  on  the  nth 
plane 

The  first  term  or.  the  right  side  of  5q  1  -resents  the  inertia  force  due  tc 
the  section  itself,  tre  second  ten  summarizes  t:e  influence  cf  the  surround¬ 


ing  water.  The  tern  vT  presents  t:  e  so-called  added  mass, 

r 


'  *-£  i  - 


assoc  iated 


witv  the  transition  of  the  r.th  section  tc  tr.e  r.  ♦  1th.  section 


F-23^5 
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The  third  term  is  a  damping  term.  The  damping  is  due  to  the  frictional 

resistance  of  the  wetted  surface,  the  generation  of  eddies,  and  the 

( 14  15 ) 

gene -at ion  of  waves  with  the  vehicle  center  as  their  origin.  '  Gen¬ 
erally  the  damping  force  is  a  quadratic  function  of  the  velocity,  tut  fol- 

( 10) 

lowing  the  method  used  for  investigation  of  ship  motions,  a  lineariza¬ 
tion  is  used.  The  effect  of  the  damping  as  a  whole  will  be  discussed  later. 

9 

The  vertical  velocity  V  and  the  vertical  acceleration  W  are  obtained 

L  D 

by  averaging  the  particle  velocities  in  the  considered  plane3  at  the  transi¬ 
tion  of  the  sections. 

Since  the  horizontal  dimensions  of  the  sections  are  small  compared  to 
the  wavelength,  the  quantitites  V  and  may  be  taken  as  the  vertical 
particle  velocity  and  acceleration  in  the  center. 


W  ~ ' 

n 


c)t> 

'z 

sr 


Acu  sinh  k(d  +  Z  -  D  ) 

-  sin  (cut  +  e) 


sinh  kd 


(4) 


z  -  Z-D 


n 


sinh  k(d  +  Z  -  D  ) 

-  cos(cut  +  e) 


a  ltd 


sinh  kd 


(5) 


z  -  n— D 


n 


where 

D  -  distance  between  the  still  water  level  and  the  average  level  of 
the  bottom  of  the  section. 

Since  tne  displacements  (Z)  are  generally  small,  compared  to  (d  -  Dp), 
E  nations  (4)  arm  (5)  reduce  to 


sinh  k(d  -  D  ) 

W  .  -  Al"> - ; — : — — —  8in(o)t  +  c) 

n  sinh  kd 


(4a) 


sinh  k(d  -  D  ) 


cos  (cut  ’  f) 


si  nit  kd 


(5a) 


r-23-5 
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For  the  vehicle  an  a  whole ,  the  linearised  equations  of  vertical  notion 


become 


•  •  re  n-1  .  n»4 

I  F  (p)  i-  £  »H.Z+M”  Z  -  I  M"  W  -  LH  W+H  Z  ♦  \  g 

n»l  n  t  n«*l  n  n-1  n  t 


where 


«t  -  EHn 


»  total  mass  of  the  vehicle 


M"  •  total  added  ness  of  the  vehicle 

H  th 

K  =  linearised  damping  coefficient  of  the  n  section 
n 

N  a  £  H 
Zt  *n 

»  external  vertical  restraining  force  of  the  vehicle 

The  total  vertical  upward  force  on  the  horizontal  planes  of  the  sec¬ 
tion  are  calculated  from  the  pres  euro  8  at  the  depths  of  these  planes.  The 
preseure  (p)  at  a  depth  z  is 


P  -  Pft-Pgz  +  PgA  -■  c~oafa  fra  co0^<0t  +  •) 


where  p  *  atmospheric  pressure 
ft 

As  the  plane b  are  "wll  compared  to  the  wavelength,  the  pressures 
on  them  are  assumed  to  be  uniform  ami  equal  to  the  pressure  in  the  center 


cf  a  plane. 


The  force  Z  F  (p)  producing  the  heaving  motion  of  the  vehicle 
u»l  a 


(Fig.  3)  is 


n»4  n»h 

“  *z  •  E.  -  pa>t  -  l  -  D  <An  -  W 

n»I  n  nal  a 


where  A  °  horizontal  cross-sectional  area  of  the  n  section 
n 


n>c4 


Using  Eqs.  (?)  aid  (8),  the  force  E  F  (~)  la  determined  as 


K-l  n 


n*4  .  n=4 

I  r ,  (P)  -  0  8  v  -  D  6  A  z  *  z  -  Vi>  C08h  11(3  *  Dn  *  z) 

n»l  n  n=l 


vhere  7  =  total  displacement  of  the  vehicle, 

x 

If  the  vehicle  is  not  restrained  (S  =*  0),  then 

r 

0  «  7t  -  Mt  s 

Since  the  displacements  (z)  are  small,  it  may  be  assumed  that 


(9) 


(10) 


cosh  k(d  -  D  +  Z)  =  cosh  k(d  -  D  ) 


(ID 


The  use  of  Eqs.  (9),  (10),  and  (ll)  in  Eq.  (6)  results  in 

n-4  .  n=4 

(M.  +  M"  )  Z  +  N  Z  +  0  g  A,  Z  =  Z  M,r  W  +  T  N  W 
t  z  '  z,  1  .zn.zn 

t  t  n=l  n  n=l  n 


g=4 

P  g  A  cos(cnt  +  e) 

cosh  kd  . 

n=l 


-  An+1)  cosh  k(d  -  Dn)  (12) 


Equation  (12)  is  in  the  form  of  the  equation  of  motion  of  a  linear 

damped  mass-spring  system.  The  mass  in  this  system,  used  as  an  analogue 

(Fig.  4),  consists  of  the  mass  of  the  vehicle  M  and  the  so-called  added 

mass  M"  ,  which  conveniently  may  be  imagined  as  an  amount  of  water  that  is 
Zt 

moving  together  with  the  vehicle.  The  damping  force  the  system  is  the 
linearized  drag  of  the  vehicle  when  moving  vertically  in  water  without  the 
disturbance  of  the  waves.  The  restoring  force  of  the  system  is  the  resort¬ 
ing  force  due  to  buoyancy,  which  is  equal  to  the  weight  of  the  displaced 
water  for  any  given  vertical  displacement  of  the  vehicle  from  its  free- 


floating  condition. 


P-2?65 
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The  excitation  of  this  analogous  system  is  represented  by  the  term 
on  the  right  side  of  Eq.  (12).  The  last  term  on  the  right  side  represents 
the  excitation  by  the  pressure  variations  cf  the  wave  on  the  horizontal 
planes  of  the  vehicle;  the  second  term  on  the  right  side  represents  the 
excitation  by  the  flow  of  the  water  around  the  body  in  the  vertical  di¬ 
rection  (drag).  The  first  term  on  the  right  side  is  difficult  to  visual¬ 
ize.  In  essence,  it  represents  an  inertia  excitation,  which  is  the  added 
mass  times  the  acceleration  of  the  water  particles  at  the  horizontal  planes. 

All  three  terns  are  periodic  forces  with  the  frequency  of  the  wave. 

It  will  be  noted  from  Eqs .  (4),  (b),  and  ( 12 )  that  these  periodic  forces 
are  not  ail  in  phase  with  the  wave.  The  magnitude  of  the  total  periodic 
force  and  its  phase  angle  relative  to  the  wave  can  be  calculated  by  a 
vectorial  summation  of  the  three  components  of  the  excitation.  Consequently, 
the  excitation  may  be  expressed  as 


n=4 

E 

n»l 


W 

n 


n=4 

+  E  N  W 
.  z  n 
n=l  n 


+ 


0  g  A  CQ8(cut  +  c) 
cosh  ltd 


n=4 

E 

n-1 


A  . )  cosh  k(d  -  D  ) 
n+1  n 


=  ?  COB  (cut  +  €  +  a)  (13) 

C  X 

where  F  =  periodic  excitation  of  the  analogous  system 
a  =  nhase  angle 

It  will  be  noted  that  the  magnitude  of  the  pressure  tern  in  Equation 
(12)  and  (13)  is  dependent  on  the  wave  frequencies  as  it  constitutes  the 
31m  of  the  pressures  on  four  surfaces.  Because  of  decreasing  pressures 
with  increasing  depths,  the  total  force  can  be  zero  for  a  particular  wave 


frequency. 
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If  no  wave  action  is  present,  the  right  side  of  Eq.  (12)  becomes  zero. 

(M.  +  M"  )  Z  +  N  Z  +  0  g  A  Z  -  0  (14; 

1  Zt  Zt  1 

Equation  (14)  presents  the  decaying  motions  in  heave  of  the  vehicle  in 
still  water  after  the  vehicle  has  been  displaced  fran  its  equilibrium 
position. 

The  motions  of  the  vehicle  in  uniform  waves  with  a  particular  height 
and  period  can  be  calculated  with  Eq.  (12)  if  the  dimensions  and  the  added 
mass  and  damping  terms  are  known. 

For  bodies  with  a  simple  geometry,  e.g.,  spheres  or  cylinders,  values 

( ip  u ) 

of  the  added  mass  can  be  calculated  by  theory. v  *  The  added  mass  of 

more  complicated  bodies  is  generally  determined  from  model  experiments  in 

still  water  by  use  of  Eq.  (14),  or— in  the  case  of  ships— i6  calculated 

from  graphs  which  are  based  on  model  experiments.  Since  the  added  mass 

M"  can  be  obtained  from  such  a  simple  experiment,  no  attempt  has  been 

made  to  derive  an  analytical  value. 

The  added-mass  terms  M"  may  possibly  be  estimated  in  proportion  to 

n 

the  areas  at  the  transition  of  the  sections  from  M"  . 

Zt 

No  damping  coefficients  were  found  in  the  available  literature  fo~ 

bodies  with  a  shape  similar  to  the  vehicle  being  considered.  By  means 

of  a  simple  test  and  u3e  of  Eq.  (14),  a  value  of  the  linearized  damping 

coefficient  N  due  to  viscous  and  wave-generating  effects  could  be  ob- 
t 

tained  for  the  vehicle  model  when  moving  in  its  natural  frequency.  This 
value,  however,  may  deviate  considerably  from  the  prototype  values,  due 
to  differences  in  Peynolds  numbers. 
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However,  since  Eq.  (12)  represents  a  damped  linear  system,  the  value 
of  damping  coefficients  is  materially  unimportant  for  the  calculation  of 
the  response  for  motions  with  a  frequency  above  twice  the  natural  frequency 
of  the  system. 

Studies  by  Keuligan  and  Carpenter  for  cylinders  and  plates  ^ ^  and 

( 17 ) 

Karleman  and  Shapiro  for  spheres  indicate  tnat  for  large  objects  the 
drag  forces  are  small  compared  with  t)"e  inertia  forces,  and  in  many  cases 
may  be  disregarded. 

Consequently,  for  frequencies  above  twice  the  natural  frequency,  E^. 
(12)  is  reduced  to 


n=4 


(Mt  +  M”  )  2  +  0  g  Al  Z  =  Z  M" 


W 


n=l  n 


n=U 


*  P  6  *  C)  t  (A  -  A  )  cosh  k(d  -  D  )  (15) 


cosh  kd 


n=l 


n  n+1' 


The  steady  state  solution  of  Eq.  (15)  is 


1  -  (— ) 

cu 

0 


1  1  034 

- - -  •  ■  Z  M"  w 

2  t  &A.  _  , 


1  n=l  n 


4  p.jg.A  ±jl  (A  .A  ,  )  co.h  k(i  -  D  ) 

cosh  kd  '  n  n+1  n 

n*l 


(16) 


where 


00 


D  g  A1 

(M  ♦  M1’  ) 
t  Zt 
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Pitch  tad  Surge 

In  the  following  analyte*  It  la  assisted  that  the  vehicle  consists  of 
a  series  of  coupled  cylinders ,  each  with  a  different  length  and  diameter. 
Kext,  It  Is  assisted  that  the  horizontal  displacement  of  each  cross  sec¬ 
tion  of  a  cylinder  is  equal  to  the  displacement  of  lta  center. 

Ualng  the  same  coordinate  system  as  is  used  for  the  analysis  of  the 
heaving  station,  and  considering  the  mass  center  of  the  vehicle,  the  equa- 


tions  of  the  surge  X  and 

pitch  | 

n-4 

n»4 

£  r  (p)  -  H,  x  - 

r  m, 

n-1  xn  x 

n-1 

0-4 

1/2 

I  c 

n-1  1 

n-4 

n-4 

I  »  (l>)  f.  •  I  * 

-  £ 

n«l  *n  n 

n-1 

n-4 

1/2 

£  C 

n-1  ( 

n  n  '  n 


(17) 


.* 


n  '  n 


^2 


n  n  n  '  n 


(18) 


where 

X  -  horizontal  displacement  of  the  mass  center 

s 

X  ■  horizontal  velocity  of  the  oass  center 

s  e 

X  -  horizontal  acceleration  of  the  mss  center 

I  -  aass  Inertia  acmrat  of  the  vehicle  around  its  mass  center 

f  «  distance  between  the  smss  center  of  a  cylinder  and  the  sms* 
D  center  of  the  vehicle 

i  a  angular  displacement  of  the  long  axis  of  the  vehicle 
• 

p  -  angular  velocity  ol'  the  long  axis  of  the  vehicle 
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M" 
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n 

M 

x 

□ 

L 

n 

R 

n 

U 

• 

u 

The 

U  = 

n 


U 
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»  angular  acceleration  of  the  long  axis  of  the  vehicle 
e  metacenter  height 
=  cylinder  number 

=  added  mass  of  a  cylinder  moving  horizontally 
=  total  horizontal  force  from  the  wave  pressure  on  a  cylinder 
=  mass  of  a  cylinder  moving  horizontally 
=  length  of  a  cylinder 
=  radius  of  a  cylinder 

*  mean  horizontal  velocity  of  the  water  mass  displaced  by  a  cylinder 

=  mean  horizontal  acceleration  of  the  water  mass  displaced  by  a 
cylinder 


mean  velocity  U  and  the  mean  acceleration  U  are  defined  aB 
n  n 


cR2 


R  x  =^R2  -  y2 
n  „  n  3 


z  =  D 


n-1 


r 

•  dz  dx  dy 


(19) 


>  2  2 

-R  x  ■  R  -  y  z  =  D 
n  n  n 


x 


D 

n 


dz  ix  dy 


(20) 


or  one  cylinder,  the  total  wave  pressure  may  be  written 
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Fx  (P) 

n 


(p  -  p  )  d2 

cl 


r‘ 


p  )  dz 
a 


x  =  1  R“  - 


where  the  wave  pressure  p  is  given  by 

p  -  pa-  p  g  I  *  0  8  A  eo‘o^i|;-^  c°s(“  1  +  £)  (22 
Evaluation  of  the  integrals  in  Eq.  (21 )  by  use  of  Eqs.  (20)  and  (22)  gives 
with  some  approxination 


r  fx  (P)  .  tor  u 

n»l  n  n=l 


n*4 

T.  F  (p)  f 


,  x 
n*l  n 


E  0  7  f  lT 
,  n  n  n 


where 


*  volume  of  the  section 


By  use  cf  Eq.  (23)  in  Eq.  ,17)  and  Eq.  (21)  in  Eq.  (i3),  and  by  linear¬ 
izing  the  damping  terras,  the  equations  of  motions  nay  be  expressed  by 


(K  -  M"  )  X  -  X  *  E  M"  f  t> 
ox  x  ,  x  n 


x  *  e  m"  r  0  -  e  f  t 

xx  ,  x  n  ,  x  r. 

t  n-^  n  r.  -  ..  n 


y;  )  -j  •  r  n 

x  n  ,  x 
- 1  n 
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n-4 


2  *, 


n-1 


n-1 


I")  +  *  L  N  f  t  +  M  g  a ,  i  +  ;■:  M"  f  X  +  £  H  f  X 

n»l  n  ^  n-1  n  n-1  n 


n— 1  ,  n—4 

£  (p  V  +  M"  )  f  U  +  -  H  f  U 

n-1  n  *n  n  n  n-1  xn  n  a 


(26) 


a.  y_ 

where  N  -  <i wiping  term  of  the  n  section  in  surge 
n 


M” 

x. 


-  total  adued  mass  In  surge 


I"  =  added-mass  inertia  mcoent 

It  will  be  noted  that  Eq.  (2^)  (pitch  equation)  represents  a  damped 
mass-spring  system.  The  excitation  of  this  system  is  coupled  with  the 
movement  in  surge. 

Following  the  theoretical  derivations  used  for  toe  heave  and  again 
disregarding  the  drag  excitation,  the  equations  of  motion  for  frequencies 
larger  than  twice  the  natural  frequency  in  pitch  are 


n— 1 


(M.  ♦  M"  )  X  +  :  M'  t  p 

X  V  n-1  *n  n 


n-1 

F  (p  V  ♦  K"  )  U 
,  n  x  n 
n-1  n 


(27) 


(I 


n— 1 
n-1 


The  solution  for  pitch  is 


n*» 4  . 

£  (p  7  ♦  M"  )  f  U 

'  a  X  QD 

n-1  n 


n-*» 

r.  m"  r 

,  x  n  n-** 
n-1  n 


1 .  h"  )  = 

t  x  '  n-1  n 


p  - 


<  \  -  2>  i  • 


n-1 

K'  f 

n-1  *n  a 
(M,  *  M"  ) 


(29) 


vt.ere 
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The  added  sms  for  a  cylinder  accelerated  in  a  direction  normal  to 
Its  axis  la  equal  to  the  vater  aaaa  displaced  by  it.^^#^  Since  for 
r.mai  t  angular  novensnts  each  section  of  the  Tehicle  la  accelerated  In  a 
direction  which  Is  approximately  in  the  direction  of  itB  axis,  it  aay 
be  as 8 used  that  I  -  I". 

For  vehicles  with  a  relatively  even  distribution  of  aass,  the  static 
sonant  of  the  added  ease  around  the  center  of  gravity  of  the  vehicle  is 


I  M"  f  «*  L  p  7  f  =>  0 

x  n  a  a 

n 


(30) 


In  this  particular  case,  the  solution  for  pitch  (Eq.  (29))  reduces  to 


^  ■ 


fn  'Jn 

_ n _ 

(Jn  -  2)  (I  ♦  I") 


£  p  7  f  U 
n  n  n 

7-2 — *T7 
(®n  -  «  )  1 


(31) 


THEORY  OF  HOVSHPTS  DPUCZD  BT  A  S1AKAY 
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In  order  to  describe  the  movente  of  the  vehicle  in  a  seevsy,  terae 
of  practical  significance  to  the  user  have  to  tot  selected.  At  present , 
descriptions  of  conplex  move  tots  are  generally  given  in  statistical  val¬ 
ues  of  the  movement j  to  conform  vith  this  practice,  this  method  of  expres¬ 
sion  has  been  used  herein. 

The  calculation  method  of  significant  values  of  responses  of  systems 

subjected  to  random  or  complex  excitations  vas  developed  by  electrical 

( 16) 

engineers  for  use  in  randan  noise  analyses'  '  and  vas  more  recently  ap¬ 
plied  to  the  calculation  of  ship  movements  in  a  seaway/^ 

The  be*.'  requirement  for  the  applicability  of  the  method  is  that 
the  response  of  the  considered  system  be  linear  to  the  excitation.  If 
the  method  is  applied  to  the  analogous  linear  system  of  the  vehicle,  as 
described  earlier  in  connection  vith  Fig.  K,  then  the  responses  of  the 
vehicle  to  a  sun  of  uniform  vave  trains  —  which  are  assumed  to  be  mov¬ 
ing  In  one  direction  --  are  equal  to  the  sum  cf  the  vehicle  responses  of 
the  individual  components . 

Since  the  responses  of  the  vehicle  in  uniform  vave  trains  can  be  de¬ 
termined  by  theory  or  by  model  experiments,  the  complex  movements  excited 
bv  vaves  vith  a  certain  spectrum  can  be  obtained  by  means  of  the  folloving 
procedure,  vhich  is  illustrated  in  Fig.  5* 

For  all  frequencies  of  the  considered  vave  spectrum,  the  value  of  the 
:nean  square  of  the  vave  solitudes  per  unit  of  vave  frequency  --  called 
spectral  denEity  of  the  waves  --  is  multiplied  by  the  corresponding  square 
of  tlm  ratio  of  response  to  eeve  amplitude,  and  the  results  are  plotted  as 


the  spectral  density  of  the  response.  In  mathematical  terms,  the  procedure 
ie  generally  written  as  follows : 

2 


6r(co )  «  Sw(<n)  ‘  |  T(to)j 


where  8  (co) 
r'  ' 


a  spectral  density  at  frequency  u>  of  the  re  spoil  8  e 


5v(a>)  *  spectral  density  at  frequency  oj  of  the  waves 


I  2 


square  of  the  ratio  of  response  to  wave  amplitude  at 
frequency  jj  of  the  wave  environment 


Once  the  spectrum  of  the  movement  in  study  has  been  obtained,  an  in¬ 
tegration  of  the  spectrum  results  in  a  single  number  E,  vhich  represents 
the  energy  of  the  movement.  If  the  spectrum  is  relatively  narrow,  the 
average  amplitudes  of  the  movement  can  be  expressed  in  terms  of  the  energy 
E,  as  shown  in  Teble  I. 


Table  1 

KESPONEE-AMPLITUIS  (A)  AND  RESPONSE -ENERGY  (E)  RELATIONS 


(Based  on  Ref.  3,  after  Ref.  1) 

Average- Response -Amplitude  Data 

Average 

aaqplitiKie  of  ail  responses 

AreB*  - 
av 

0.886 

VT 

Average 

ssqpiitude  of  the  1/3  highest 

resp 

*1/3  " 

1.416 

V* 

Average 

amplitude  of  the  1/10  highest 

a5^bp  - 

l.8oo 

vr 

Greatest-Re sponae-Amplitude  Probability  Data 

Number  of  Nine  times  out  of  ten,  tne  highest  amplitude  of  the 
cycles  (S)  response  out  of  a  series  of  N  cycles  is  between 


20 

1.40 

Ik. 

and 

2.44 

7B 

100 

1.38 

"V  E 

2.7^ 

Ve 

‘oOQ 

2.26 

VI 

3.03 

V'E 

1000 

2.41 

V  E 

3.14 

",  £ 

Note;  The  total  deviations  of  the  response  from  minimum  to  aiajciuum 
are  twice  the  values  of  the  amplitudes  in  the  table 
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TEST  FACILITIS6  ASP  FROCZDURBS 

The  experimental  investigations  were  performed  in  a  water-wave  tank 
100  ft  long,  k  ft  high,  and  2  ft  vide. 

At  one  end  of  the  tank,  waves  could  be  generated  by  a  bulkhead-type 
veve  generator.  A  gravel  beach  vith  a  1-10  slope  at  the  other  end  of  the 
tank  vas  installed  to  absorb  the  generated  waves  and  to  reduce  the  genera¬ 
tion  of  stanulng  vaves. 

The  vehicle  models  used  in  the  experimental  investigation  were  wooden 
models,  17  ~  in,  long -and  lj-  in.  in  diameter,  representing  a  1-120  scale 
model  of  a  million-pound  solid-propallaat  rocket  vehicle,  175  ft  long  and 
12  £  ft  in  diameter. 

The  movements  of  one  of  the  models  were  sensed  by  use  of  tvo  displace¬ 
ment  gages  (Pig.  6).  Each  gage  consisted  of  a  thin  wire  soldered  to  tvo 
very  thin  brass  plates,  one  located  vertically,  the  other  horizontally. 

Strain  gages  were  mounted  on  these  plates.  Movements  of  the  ends  of  the 
wires  (arms)  of  the  displacement  gages  induced  strains  in  the  strain  gages, 
which  were  amplified  by  carrier  amplifiers  and  recorded  simultaneously  with 
the  output  of  a  resistance-type  water-level  recorder  in  a  high-speed 
oscillograph. 

It  vas  considered  that  the  displacement  gages  did  not  restrain  the 

* 

Eodel  in  its  orbital  movements,  since  the  forces  necessary  to  displace  the 
end  of  the  gages  were  insignificant  (approximately  50  mg  par  inch  deviation). 

t 

The  movements  of  the  model  during  experimental  test  runs  were  recorded 
from  the  moment  the  first  ’rare  reached  the  model  until  obvious  irregulari¬ 
ties  due  to  standing  waves  were  generated  in  the  tank. 


To  obtain  information  about  the  characteristics  of  the  vehicle  in 
free  hears  and  la  free  pitch,  the  model  vas  displaced  vertically  and  hori¬ 
zontally  and  suddenly  released  after  each  displacement.  The  movements  were 
recorded  until  they  had  decayed  to  insignificant  levels. 

To  check  on  the  results  obtained  vith  the  displacement  gages,  photo¬ 
graphic  observations  vere  made  vith  a  second  model.  This  model  was 
equipped  vith  tvo  small  light  bulbs,  one  at  the  head  and  one  at  the  tail. 
Power  to  these  bulbs  vas  supplied  by  means  of  tvo  vires  0.005  in.  diameter. 
(These  vires  were  considered  not  to  restrain  the  motions  of  the  model.)  By 
means  of  a  mechanically  operated  switch,  the  light  bulbs  flashed  on  for 
0.02  sec  every  0.15  sec.  By  means  of  simple  electrical  circuits,  the  flashes 
and  the  time  of  exposure  of  the  photographic  film  in  the  camera  vere  recorded 
which  vas  located  near  the  model.  The  scale  on  the  photographs  vas  obtained 
by  use  of  white  vine  grids  placed  at  both  sides  of  the  tank  (Pig.  1) . 
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III.  RESULTS 

UNIFORM  WAVS  TRAINS 
Experimental  Results 

Figure b  7  and  3  present  all  experimental  results  as  functions  of  the 
wave  frequency.  (Theoretical  results,  which  will  be  treated  later  in  this 
report,  are  plotted  for  comparison.)  Figure  7  shows  the  results  of  all 
lata  analyses  for  heave.  The  data  are  presented  as  plots  of  the  ratios  of 
double  heave  amplitude  to  wave  height  (S^/H) ,  which  equals  the  ratio  of 
heave  amplitude  to  wave  amplitude  (A^/A)  as  a  function  of  the  wave  frequency 
and  period.  It  will  be  noted  that  the  ratios  increase  with  a  decrease  in 
frequency. 

Figure  6  shows  the  results  of  all  data  analyses  of  the  maximum  angular 
deviations  (pitch)  from  the  vertical.  The  data  sue  presented  as  plots  of 
the  ratio  of  the  maximum  angular  deviation  in  minutes  of  arc  to  the  wave 
height  in  tenths  of  an  inch  as  a  function  of  the  wave  frequency.  The  cor¬ 
responding  prototype  frequencies  and  periods  are  Indicated  at  the  top  of 
the  graph.  The  prototype  ratios  are  expressed  as  the  maximum  angular  devia¬ 
tion  in  minutes  of  arc  to  the  wave  height  per  foot. 

Figure  9  presents  the  results  of  a  photographic  observation.  The  shut¬ 
ter  of  the  camera  was  opened  for  approximately  one  second.  The  time  of  the 
opening  vac  recorded  by  the  oaciHograpu,  together  with  the  time  of  the 
light  flashes  on  tJie  need  and  tail  of  the  model.  The  outpvot  of  the  wave 
recorder  was  also  recorded.  The  positions  of  the  head  and  tail  of  tne 
model  cen  be  identified  in  the  photograph,  which  is  reproduced  in  Fig.  9. 

In  order  to  snow  the  positive  of  the  vehicle  in  more  detail,  the  positions 
■f  the  head  aau  tail  are  enlarged  ten  times,  ir»d  lines  are  drawn  represent¬ 
ing  the  movementc  of  the  long  axis,  enlarged  ten  times. 
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Xo  measurements  wn  made  of  svay,  y*w,  and  roll.  Visual  observations 
of  tbs  aoT—euts  of  tbs  vsbiols  in  uniform  vavss  in  tbs  wave  tank  indicated 
that  tbsss  motions  wars  negligible. 

Within  tbs  range  of  tbs  wave  heights  used  in  tbs  experimental  inves¬ 
tigations  equivalent  to  3-  to  7 -ft -high  vavss  for  tbs  prototype,  tbs  re¬ 
sponses  of  tbs  model  in  heave  wad  pitch  vers  significantly  linear  vith  the 
vavs  height.  Tor  example,  it  is  estimated  that  the  experimentally  obtained 
ratio  of  heave  to  vavs  height  for  a  wave  height  of  6  ft  in  prototype  is  ap¬ 
proximately  10  per  cent  smaller  than  the  ratio  for  a  vavs  height  of  3  ft 
in  the  prototype. 

Theoretical  Results 

The  added  mesa  for  '  uses  calculations  was  determined  by  use  of  the 
results  of  a  modal  test  in  still  vater.  In  this  test,  the  model  was  dis¬ 
placed  vertically  end  released.  The  following  decaying  motions  were  meas¬ 
ured  and,  by  use  of  Iq.  ( 15) ,  the  added  mass  calculated. 

The  results  of  the  calculations  for  tbs  bearing  mot  ion  vers  obtained 
by  use  of  Xq.  (l6e).  The  experimental  results  presented  in  Tig.  7  pertain 
to  vehicle-model  movements  in  a  vater  depth  of  2.68  ft  (prototype  320  ft  ). 
Since  the  particle  motions  of  waves  with  a  length  larger  than  twice  the 
vater  depth  are  affected  by  the  bottom,  the  response  in  deep  vater  (d  >  l/2X) 
will  deviate  from  the  response  at  a  depth  of  2.68  ft  (prototype  320  ft). 

The  results  of  tbs  theoretical  stuty  of  the  heave  are  in  good  agree- 
merit  vlth  the  results  of  the  tests  made  in  the  wave  tank.  The  maximum  de¬ 
viation  of  the  theoretical  values  from  the  average  values  obtained  by  ex¬ 
periments  is  approximately  20  per  cent. 
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The  results  of  the  calculation  for  pitch  for  she  model  were  »tde  by 
use  of  Eq.  (29)  and  platted,  in  Pig.  3.  The  oalcvlahiocs  were  made  by  using 
tue  experimentally  determined  natural  freq  ency  in  pitch.  Again  the  theore¬ 
tical  results  are  in  fair  agreement  with  he  experimental  results. 

Seaway  -  Theoretical  Results 

As  an  illustration  of  the  application  of  spectrum  calculations,  the 
results  for  heave  and  pitch  of  the  prototype  vehicle  acted  upon  by  nonuni- 
form  waves  in  deep  water  are  pro  seated.  :.n  Table  II.  A  fully  arisen  sea  in¬ 
duced  by  wind  speeds  of  20  kn  was  selected  as  a  typical  seaway.  Thia  wave 
environment  represents  conditions  for  a  Sea  State  k  (Hydrographic  Office 
Scale)  without  swells. 

The  calculations  were  made  by  use  of  the  results  of  the  theoretical 
study  of  heave  of  tie  vehicle  in  uniform  deep-water  waves  (Pig.  7)  an*1 
the  results  of  tie  experimental  study  for  pitch  (Fig.  8).  Figure  10  pre¬ 
sents  the  calculations  according  to  the  pattern  presented  in  Pig.  6  for 
the  vehicle  movement  in  heave. 


Table  a 


smvua  or  calcuxatkm  for  hkavi  akb  pitch 

II  A  FULL?  AftJSfif  SKA 


i 

1 

1 

r  ■ 

Fully  Arisen  Sea  with 
Wind  Speed  of  20  kn 
(Sea  State  k) 

, 

Characteristic 

Statistical 

Value 

Heave 

(double 

aoplltude, 

ft) 

Pitch 

(aaxlmon  devia¬ 
tion  frcm  the 
vertical,  min 
of  arc) 

Average 

j  0.8 

50 

Average  1/Z 

highest  ! 

( significant) 

1.3 

j 

t 

1 

80 

» 

Average  l/lO 
highest 

i 

■ 

1.7 

j 

100 

1 

In  seriea^of  1000  | 

cycles,*1  9  times  1 
out  of  10  the  ' 

highest  values  • 
Is  between 

2.2  -  2.9 

i 

1 

1 

130-110 

1 

*For  a  fully  arisen  im  with  a  vind  speed  of  20  fcn,  the 
respective  time  for  1000  cycles  is  appmtlmtoly  2  {_  hr. 


-jHu,*--.. 
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IV.  DISCUSSION 

It  appears  from  Figs.  7  and  8  that  the  experimental  results  are  in 
fairly  good  agreement  vith  the  theoretical  results.  Discrepancies  may¬ 
be  explained  by  the  approximations  made  in  the  derivations  of  the  formu¬ 
las  used  and  in  the  accuracy  of  the  measurements.  A  higher  accuracy  will 

be  obtained  by  considering  the  movements  of  the  two-dimensional  (horizontal) 

.  .  f oq  ) 

flow  around  cylindrical  sections  of  the  vehicle  with  a  very  snail  hcight(dz).  '  ' 

Since  the  particle  motions  in  shallow  water  (prototype  wave  periods 
larger  than  11  sec)  are  affected  by  the  bottom,  the  orbits  are  ellipti¬ 
cal  rather  than  circular;  consequently,  the  movements  of  the  vehicle  in 
shallow  water  deviate  from  the  movement  in  deep  water,  which  is  illustrated 
by  the  heaving  motion  in  Fig.  7. 

Since  the  particle  motions  near  the  surface  are  larger  than  those  at 
a  great  depth  and  the  excitations  are  related  to  the  volumes  and  their 
accompnaying  added  masses  (Eqs.  (13)  and  (29)),  a  reduction  of  the  dimen¬ 
sions  near  the  surface  will  tend  to  reduce  the  movements. 

Smaller  dimensions  at  the  water  line  will  also  generally  decrease 
the  frequency  of  free  heave  and  pitch,  and  the  chance  for  resonance  with 
low-frequency  components  of  the  wave  spectra  will  be  reduced.  Since  long 
and  slender  vehicles,  as  investigated,  have  a  low  damping,  waves  with  fre- 
quen  es  close  to  the  frequency  of  free  heave  may  induce  excessive  movements. 

Consequently,  swells,  which  are  waves  generated  at  s.  considerable  dis¬ 
tance  from  the  point  of  observation,  may  induce  the  raa^or  heaving  motion 
for  large  vehicles. 


Such  movements 
re  ac^ca; itched  by  1 
face  of  the  vehicle, 


can  be  decreased  by  increasing  the  damping.  This  can 
cstalling  horizontal  plates  or  fins  on  ti.“  outer  aur- 
preferatiy  close  to  the  tail  end.  In  addition  to  an 
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inortM*  In  tapping,  such  plates  also  vill  Increase  tha  added  mas  a,  and 
consequently  a  favorable  lover  frequency  In  free  heave  vill  be  obtained. 

The  application  of  the  derived  formulas  for  heave  and  pitch  is  per¬ 
mitted  for  frequencies  larger  than  tvioe  the  frequencies  of  free  heave 
and  pitch.  Thus,  in  ths  oase  of  the  studied  vehicle,  spectra  calculations 
aa  presented  in  Fig.  10  suet  be  limited  to  seaways  vhich  contain  no  fre¬ 
quency  cosaponeirtg  lover  then  tvice  tha  frequency  of  free  heave  (periods 
larger  thar.  3  sec). 

In  the  spectra  calculations  it  is  assumed  that  all  vaves  are  propa¬ 
gating  in  one  direction.  In  reality,  however,  vaves  are  generated  in  dif¬ 
ferent  directions,  hut  generally  a  very  small  part  of  the  vmve  energy  is 
in  eaves  vhich  make  a  significant  angle  (larger  then  approximately  15  de¬ 
grees)  from  the  main  direction/ ^  Because  of  the  small  horizontal  di¬ 
mensions  of  the  vehicle,  the  heave  vill  not  be  affected  by  directional 
effect  a;  tho  pitch,  however,  is  dependent  on  the  directional  component  c 
of  the  spectrum. 

It  seems  possible  to  determine  the  novenents  of  a  vehicle  in  ocean 
vaves  by  calculating  the  responses  at  the  different  directions  of  the 
spectrum  and  a  vectorial  addition  of  the  cocgjonenta  thus  obtained.  Since 
the  spread  of  the  main  veve  direction  is  relatively  small,  the  spread  of 
the  responses  is  small  and  its  effect  can  be  disregarded. 


V.  COHCUJSICBS 


The  ncrvenents  In  heave  and  pitch  of  a  large  rocket  vehicle  drifting 
in  a  vertical  attitude  in  unifona  vavee  can  be  calculated  by  naans  of  the 
t}»ory  preaented  in  tbit  paper  if  the  periods  of  the  vavea  are  shorter 
than  half  the  periods  of  free  heave  and  pitch.  The  result*  of  such  cal¬ 
culations  for  unifora  waves  can  be  used  for  the  calculation  of  the  fre¬ 
quency  of  occurrence  of  certain  levels  of  the  heave  and  pitch  response  of 
the  vehicle  to  a  seaway  vlth  a  given  wave  spectrum* 
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Fig.  i  —  One-to-120  scale  models  of  the  rocket  vehicle 
in  the  wave  tank  (  observation  grid  width  is  2  in.) 
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Fig.  3  —  Approximation  of  the  model  dimensions 

for  the  analyses 


p 


\ 


\ 


/ 

/ 


St 

CO  H 


cc 

Q_ 


>- 

cr 

<r 


co 

co  z: 
<X  c  o 


c_> 

X 

LU 


> 


c 


II 

<t 

cr 

Q_ 

+ 

■rx 


cp 

Ll_ 


Mechanical  analogy  of  the  vehicle  in  heaving  motion 
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Fig 


b  -- Calculation  of  a  vehicle  -  response 
spectrum  Fjm  a  giver  wave  spectrum 


TOP  OF  WAVE  TANK 


P-2365 

35 


Z 

o 

5 

cr 

CO 

< 


UJ 

m 

D 

h- 

5 

Z) 


_I 

UJ 

Q 

O 

'i 


z  y 

<  >— 
cr  ct: 

t_  UJ 

CD  > 


O 
cr  cr 

H  O 
c n  x 


to 

<D 

cn 

o 

cn 


c 

<U 

E 

Q> 

O 

CL 

CO 

*0 

(1) 

_c 


o 

* 

a> 

> 

o 

o 

E 

0) 

_c 

o 

CO 


CD 

& 

Li. 


DOUBLE  HEAVE  AMPLITUDE/WAVE  HEIGHT  (Hh/H) 
HEAVE  AMPLITUDE/ WAVE  AMPLITUDE  (Ah/A) 


P-P365 

36 


r - • - 1 - * - 1 - - - 1 - - - » - > - r - ^  U) 

0.2  0.4  0.6  0.8  1.0  1.2 

CIRCULAR  FREQUENCY  (RAD/SEC) 

—  T  -  — T  - - - - I - 1 - T - 1 - - ? - T 

20  15  10  9  8  7  6  5 

PERIOD  (SECONDS) 


50  20  10  05 

PERIOD  (SECONDS) 

Fig  7  —  Ratio  of  heave  to  wave  height  as  a  function  of  the 
wave  frequency  for  the  model  in  the  wove  tank  and 
the  prototype  in  320- ft  and  in  deep  water 


prototype 

SCALES 


MAXIMUM  DE  VIATION  FROM  VERTICAL/WAVE  AMPLITUDE 
(MINUTES  OF  ARC/FT  FOR  PROTOTYPE) 


P-2365 

37 


50 


40 


30 


20 


10 


0 


PROTOTYPE 

SCALES 


F- 

X 

o 


25 


UJ  Ld 

>  Q  on 

<  O 
*  5 

\ 

_i  cr 
<1  o 

U  U_ 

or  5  15 

UJ  Z 

>  - 


5 

O 

a: 


o 

cr 

< 


^  O 

t- 

<  ‘St 
UJ 


> 

UJ 

Q 

5 

X 

5 


z 


10 


5 


X 

< 


5 


0  L 

0 


MODEL 

SCALES 


0  2  '  04  0^6  08  To  '  r  2  ^ 

CIRCULAR  FREQUENCY  (RAD/SEC) 

20  15  10  9  8  7  6  5 

PERIOD  (SECONDS) 

°  MEASUREMENT  MADE  WITH  DISPLACEMENT  GAGES 
*  MEASUREMENT  BY  PHOTOGRAPHIC  MEANS 
v  CALCULATED 


o 


DEEP  WATER  CONDITION 

F - 

IN  WAVE  TANK 


DRAG  CHARACTERISTICS  OF  VEHICLE 
ARE  OF  MINOR  IMPORTANCE 


2  4  6  8  10  12  14 

CIRCULAR  FREQUENCY  (RAD/SEC) 

30  20  10  C 5 

PERIOD  (SECONDS) 


u’ 


Fig  8 -- Ratio  of  maximum  deviation  from  the  vertical  to 
wave  height  as  a  function  of  the  wave  frequency 
for  the  model  in  the  wave  tank  and  the 
prototype  m  a  320- ft  water  depth 
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